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The notion of an atom-light quantum interface has been developed in 
the past decade, to a large extent due to demands within the new field 
1 of quantum information processing and communication. A promising 

type of such interface using large atomic ensembles has emerged in the 
past several years. In this article we review this area of research with 
a special emphasis on deterministic high fidelity quantum information 
protocols. Two recent experiments, entanglement of distant atomic 
objects and quantum memory for light are described in detail. 
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1 Introduction 

-i— > ■ 

Interaction of light with atoms has been always one of the most exciting ar- 
eas in AMO physics. A complete quantum-mechanical approach to it allows 
for considering such aspects as quantum state transfer between atoms and 
light and generation of entangled states of light and atoms. Entangled or 
non-separable states form the basis for Quantum Information Processing and 
Communication where information is encoded, processed and transferred in 
quantum states of light and matter (for a collection of review articles on QIPC 



see, for example, Zoller et al. ( 20051 )). Processing information encoded in quan- 



tum states provides parallelism and unconditional security, two most promis- 
ing properties within this field. State transfer between matter where quantum 
information is processed and stored, and light, the prime long distance carrier 
of information, is crucial for quantum networks and other applications. 

A well-designed interaction between two or more quantum systems is an effi- 
cient tool for creation of a desired joint entangled state of the systems. Such 
interaction alone, or combined with a projective measurement on one of the 
systems can lead to the creation of a non-trivial quantum state of the remain- 
ing system(s). In this review we describe recent progress in the development 
of a tool box for deterministic quantum state generation of light and atomic 
ensembles via a well known dipole interaction. This tool box forms the basis 



1 



for a quantum atom-light interface capable of generating entangled states of 
atoms and light on demand, performing efficient quantum state exchange be- 
tween light and atoms - quantum memory for light, and long distance transfer 
of a quantum state of atoms via teleportation and other quantum transfer 
protocols. 



A key to successful quantum state engineering is a favorable balance be- 
tween efficient unitary interaction between the systems in question (light and 
atoms within the context of this article) and decoherence caused by interaction 
with the environment. Historically the first setting considered most promising 
for the light-atoms interaction was Cavity Quantum Electrodynamics (Cav- 
ity QED) where the Jaynes-Cummings type of interaction could be realized. 
Drive n by the interest in single photon - single atom interaction ( Cirac et all 



19971). this field of research has achieved spe ctacular success es in both opti- 



cal ( VlcKeever et al. , 2004 ; Kuhn et al. , 20021 ) and microwave ( Raimond et al 



2001) domains. Through the dramatic improvements in cavity quality factors 



and with the use of cold atoms, coupling of single atoms to single photon cav- 
ity modes has been demonstrated. However, formidable technical challenges of 
the strong coupling cavity QED regime promoted search for alternative routes 
to atom-light quantum interfaces. 



One of the most successful of these emerged in the second half of the 90s 
with the recognition of the fact that the use of collective quantum states of 
large atomic ensembles instead of single atoms can provide efficient "strong" 
coupling to light outside of the cavity QED regime. An intuitive explanation 
goes as follows. In order to be efficiently coupled to a single atom, light has 
to be focused and shaped to match the atomic dipole pattern. Since this 
is very difficult in free space, a cavity enhancing interaction of light with 
a specific mode is necessary. On the other hand a collection of atoms with 
transverse dimensions large compared to the wavelength of light collectively 
couples very efficiently to any spatial light mode which matches its shape. In 
classical electrodynamics this is manifested by large absorption and dispersion 
that characterizes an atomic ensemble interacting with nearly resonant light. 
Luckily, the same scaling persists for some collective quantum properties of 
light and atoms. 



The first approach to coupling quantum features of light to an atomic ensemble 
was based on a straightforward assumption that if the input light is completely 
absorbed by atoms, its quantum f eatures should be to a certain extent trans- 



ferred onto atoms. The proposal (|Kuzmich et all Il997l ) considered a V-type 



atomic level scheme where excitation in the two arms of the V-transition is 
carried out by quantum correlated - entangled - light modes. Absorption of 
these modes leads to the creation of entanglement in the excited state of atoms. 
If excitation, for example, is performed with squeezed light, the excited state 
of atoms becomes spin squeezed. In case of cw excitation spontaneous emis- 
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sion was shown to limit the degree of spin squeezing of atoms to 50%. The 
experimental implementation of this proposal has been carried out for cold 
cesium atoms excited by squeezed light resulting in the fi rst demonstration of 
a macroscopic (10 6 atoms) ensemble of entangled atoms (Iflald et al. . 20051 ). 



Efficient coupling at the quantum level also requires that no other modes of 
light but the desired input one couple to the atoms. In classical physics other 
modes are of no interest provided they do not contain any light. In quantum 
physics any empty mode is in a vacuum state. Coupling to such modes leads 
to spontaneous emission of atoms, the major source of decoherence which a 
quantum state engineer has to fight. In the Cavity QED setting interaction 
with vacuum modes of the field is made small compared to coupling to the 
relevant input mode by the strong coupling provided by the cavity. How can 
this type of decoherence be suppressed for a large atomic ensemble in the 
absence of any cavity? The solution is to use for coupling with light two 
substates of the atomic ground state rather than a ground and an excited 
state as in the standard Jaynes-Cummings approach. The advantages of using 
ground state levels are multiple. Their spontaneous decay is negligible leading 
to long coherence times. Their energy spacing is in the radio-frequency or 
microwave domain which means that with the size of the ensemble smaller than 
the radio- or microwave- wavelength the position of each atom is irrelevant 
and, hence, collective coupling insensitive to atomic motion can be achieved. 



In order to circumvent limitations due to spon taneous emission, the idea of 
complete absorption has been taken further by Kozhekin et al. ( 2000f ) where 
a driven Raman transition involving a weak quantum mode in one arm has 
been shown to be capable of faithful transfer of the light quantum state 
onto ground state collective atomic coherence. In a parallel development, a 
celebrated electro-magnetically induced transparency process utilizing care- 
fully timed Raman puls e s has been proposed for quantum storage of light 
( Fleischhauer and Lukin . 2002[) . The first exper iment testing these ideas has 
been recently carried out ((Eisaman et al. , 2005h . 



In all of the above approaches to light-atoms quantum interface the tool 
box has only included light-atoms interactions. A significant next step al- 
lowing implementation of atomic entanglement and quantum memory has 
been made by adding a quantum measurement and feedback to the pic- 
ture. The two experiments which are central to this article, the determin- 
istic entanglement of distant atomic objects and the deterministic quantum 
memory for light, have been carried out following essentially the same sce- 
nario. A pulse of light interacts with atoms via a quantum non-demolition 
(QND) Hamiltonian, a projective homodyning measurement is performed on 
the transmitted light, and the feedback conditioned on the result of the mea- 
surement is applied to atoms (Fig. |T]and|2J). Early proposals for using the 
QND interaction for atomic quantum state engineering have been published 
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by Sanders and Milburn ( 19891 ) . Cornev and Milburn ( 1998I ). and Wiseman 
( 1998f ). In particular, spin squeezing in an atomic ens emble generated by a 



QND measurement has been studied theoreticall y in ( Kuzmich et al. . ,1998; 



strated in (Kuzmich et al., 200 



Takahashi et al.l.ll99 9: Madsen and Molmerl.12 004) and e xperimentally demon 



Geremia et all |2004J). Quantum feedback 



has been explored in a series of papers by Wiseman and co-authors. In par- 
tic ular quantum feedbac k in relation to spin squeezing has been analyzed 



in (|Thomsen et all 120021) 



To complete this brief introduction into atom-light quantum interface research, 
we wish to draw a line between the deterministic interface reviewed in this 
article and probabilistic schemes. The latter are based on generation of entan- 
glement between light and atoms conditioned on a random event of detection 
of a photon in a certain mode. Spectacular experimental progress along these 



lines has bee n achieved for schemes involvin g a single ion ([Blinov et all 12 004; 



Polzikl.l2004[), single atom (IVolz et all (12005^. a nd atomic ensembles (lK uzmic h et al.. 
2003 [lvan der Wal et all 120031: IChaneliere et all 120051: IChou et al! 1200,4 lEisaman et al 
2003). 



The paper is organized as follows. In section 2 we introduce quantum variables 
for light - Stokes operators, and for atoms - collective spin operators. An effec- 
tive Hamiltonian and equations of motion for the atomic and light variables 
are then derived. We then expand the theory to two atomic ensembles with 
an addition of a bias magnetic field. 

In section 3 we formulate the equations of motion in the language of canon- 
ical variables for light and atoms. Two main quantum information protocols 
are then described theoretically: entanglement of two atomic ensembles and 
quantum memory for light. In both case the treatment includes a simple de- 
coherence model. The quantum memory protocol which we dubbed "the di- 
rect mapping protocol" represents an alternative to teleportation method of 
a high-fidelity quantum state transfer from one system to another. Unlike 
teleportation it involves only two objects: one carrying the input state and 
another the target on which the input state is transferred. 

Section 4 describes experimental methods and the key elements of the setup 
including atomic cells, magnetic fields and detection of light. Application of the 
magneto-optical resonance method for characterization of the collective atomic 
spin states is outlined. Magnetic feedback used for the spin manipulation is 
discussed. 



Section 5 contains the main experimental results on atomic entanglement and 
quantum memory for light. The determination of the benchmark atomic quan- 
tum noise level - the projection noise - is described in detail. The effects leading 
to decoherence of generated atomic entangled states are summarized. 
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In the appendices we discuss various technical aspects, such as the effect of 
atomic motion on quantum state generation and decoherence, the influence of 
polarization of light, and of technical noise of the probe laser. 

We conclude with a brief summary and outlook. 



2 Atom-Light Interaction 

In this section we introduce the quantum variables for light and atoms and 
describe the off-resonant dipole interaction between the ground state 6Si/ 2 
and the excited state 6P 3 / 2 in cesium. We will use spin operators for atoms 
and Stokes operators for light as a convenient way do describe the interaction, 
and we will present an effective Hamiltonian which describes the dynamics 
of the ground state spin and the light. For simplicity the Hamiltonian will 
be specific to the cesium ground state, although, the same procedure can be 
applied for any other atom with a magnetically non-degenerate ground state. 
With this as a starting point we derive equations of motion for the light and 
atomic operators of a single atomic ensemble. 

A significant next step for quantum information protocols discussed in this 
article is the transition from a single atomic ensemble to two oppositely ori- 
ented ensembles accompanied by the addition of a constant magnetic field 
(Fig. . If this field is oriented along the mean atomic spin direction it allows 
for achieving entanglement between two ensembles with a single pulse of light 
instead of two pulses required in the absence of the field. An even more impor- 
tant advantage is brought about by the possibility to make measurements at 
a rather high Zeeman frequency, thus achieving quantum limits of sensitivity 
with macroscopic numbers of atoms via spectral filtering of classical noise. 

2.1 ATOMIC SPIN OPERATORS 

The ground states of cesium are characterized by its outermost electron which 
is in the 6S1/2 state, i.e. the orbital angular momentum L is zero. The electron 
spin S and thus the total electronic angular momentum J has quantum number 
S = J = 1/2. The nuclear spin I of cesium-133 has / = 7/2, and the coupling 
between the nucleus and the electron gives rise to the total angular momentum 
F = I + J with quantum numbers F = 3 and F = 4. 

It is indeed the total angular momentum F which interests us in this work 
since F and the magnetic quantum numbers mp define the energy levels of the 
ground states in the limit of low magnetic field discussed here. Furthermore, 
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Figure 1. Schematic layout illustrating main components of the experiments. Two 
clouds of Cs atoms contained in paraffin coated glass cells are placed inside the 
magnetic shields. A light pulse consisting of orthogonally polarized strong classical 
and weak quantum fields passes through the cells and is detected by the detector on 
the right. Electronic feedback is applied to the rf magnetic coils surrounding cells 
to rotate the collective atomic spin. The bias magnetic field is applied to atoms in 
both cells. Atoms in the two cells are optically pumped as shown in the figure. 

we restrict ourselves to one hyperfine level, F = 4, which is possible experi- 
mentally since the hyperfine splitting = 9.1926GHz is large compared to 
typical resolutions of our laser systems. We choose to denote the total angular 
momentum of a single atom by j and for a collection of atoms (in the F = 4 
state) we denote the collective total angular momentum by J, i.e. 



where N is the number of atoms in the F = 4 state and j(*> is the total angular 
momentum of the z'th atom. The reason for using J and not F is conventional, 
and indeed, we wish to think about our spins more abstractly than just the 
properties of some atoms. Many of the results should be applicable in a broader 
sense than to a collection of cesium atoms. 

In our experiments the number of atoms N is of order 10 12 and we will almost 
always aim at having all atoms polarized along one direction which we denote 
as the x-axis. With the x-axis as quantization axis we have mp = 4 for all 
atoms to a high degree of accuracy, and the collective spin J x will really be 



N 




i=l 
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a macroscopic entity. With this experimental choice, we may treat the x- 
component of the collective spin as a classical c-number, i.e. we replace the 
operator J x by the number J x . The transverse spin components J y and J z 
maintain their quantum nature. They will typically have zero or a small mean 
value. The quantum fluctuations are governed by the commutation relation 
and the Heisenberg uncertainty relation (with h — 1) 



Jyi Jz 



Var(J„).Var(J,)>5 



(2) 
(3) 



With 10 12 atoms the quantum uncertainty of the angle of the collective spin 
direction is of order 10~ 6 . 



2.2 POLARIZATION STATES OF LIGHT 



All our experiments involve narrow-band light interacting with atomic spin 
states, and it turns out that the polarization states of the light form a conve- 
nient language to describe the light degrees of freedom. 

Consider a pulse of light, or a collection of photons, propagating in the z- 
direction. The polarization state is well described by the Stokes operators 

S x = - (n ph (x) - h ph (y)) = - (a ] x a x - a ] y a y ^ , 

S y = ^ (n ph (+45°) - n ph (-45°)) = - {a\a y + a\a x ) , (4) 

&* = \ ("pt(°+) - V0-)) = Yi ~ h l hx ) ' 

where n p h(x) is the number of photons in the pulse with x-polarization, and 
so on. The Stokes operators are dimensionless as written here, they count 
photons. At our convenience we will later break these up into time or spatial 
slices. 

We make the following assumption central to the experiments described below. 
We assume that light consists of a strong component linearly polarized along 
the x-direction and a much weaker component polarized in the y-plane (Fig. [J). 
This means that we can treat the x-mode operators a x , a x , and hence S x — > 
S x as a c-number. Note, this is very similar to the approximation of a well 
polarized sample of spins in the previous section. Specifically we find (with 
a x = a x = A x being real) that S y = A x ■ (a y + a y )/2 and S z — A x • (a y — a y )/2i. 
We see that in our approximation the quantum properties of S y and S z are 
solely encoded in the ^/-polarized part of the light. 
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It can be shown that the Stokes vector satisfies angular commutation relations 



S y , S z — iS x 



Var(S y ) ■ Vai(S z 



C2 

> £k 

~ 4 



(pulse of light). 



(5) 
(6) 



2.3 OFF-RESONANT COUPLING 



We consider the special case with a propagating beam of light coupled off- 
resonant ly to the 65'i/2 i f=4 - ► 6P3/ 2) F'=3,4 i 5 transitions in cesium. Then ab- 
sorption effects can be neglected and all the dynamics are of dispersive nature. 
Furthermore, the optically excited states can be adiabatically eliminated and 
an effective Hamiltonian describ ing the light interacting with only ground s tate 
degrees of freedom is obtained ( Julsgaard, 20031 : Hammerer et al. . 2005a ) 1 . 



cff 
int 



hey A 2 
8AA2tt Jo 



a ■ 4>(z, t) + oi • S z (z, t)j z (z, t) 



+ a 2 



{z, t)f z (z, t) -3-(z, t)f + (z, t) - S + (z, t)f_ (z,t) pAdz. 



(7) 



Here we have assumed a one-dimensional theory for the light which is sufficient 
for a beam cross section A that is much larger than the squared wavelength 
A 2 . The small letter spin operators j(z, t) are dimensionless and refer to single 
atoms at position z at time t. The Stokes operators S(z, t) are taken to count 
number of photons per unit length at position z and time t. The integration 
then runs over the entire sample of length L with atomic density p. In the front 
factor 7 is the natural FWHM line width of the optical transition 6S1/2 — > 
6P3/2 and A is the detuning from the F = 4 to F' = 5 transition with red 
being positive. 

As for the operators, 4>(z, t) is the photon flux per length, S+ — S x + iS y = 



and 5*_ 



St. iS v 



;t 



a + are raising and lowering operators 



converting <7 + -photons into a_-photons or vice versa, j± — j x ± ij y are the 
usual raising and lowering operators for the spin. 

The parameters a , ai, and a 2 are for the F = 4 ground state in cesium given 



^ (jjulseaardl 2003t ) contains a factor of two error which has been corrected in 
Eq. ©. 
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by 



a = T I ~, X 7T + 



4Vl-A 35 /A I-A45/A 

1/35 21 



a 1 



a 2 



120 V 1 - A35/A 1 - A45/A 

1/5 21 
240 ll- A35/A " 1- A45/A 




where the limit is calculated for very large values of the detuning. The de- 
tunings A 35 /27r = 452.2 MHz and A 45 /27r = 251.0 MHz are given by the 
splitting in the excited state. Let us comment on the different terms in the 
Hamiltonian (J7J). The first term containing ao just gives a Stark shift to all 
atoms independent of the internal state but proportional to the photon density 
<fi(z, t). The second term containing a± rotates the Stokes vector S and the spin 
vector J around the z-axis, known as Faraday rotation - we look more closely 
into this below. The last terms proportional to a 2 are higher order couplings 
between the light and the atoms and since a 2 is small for a sufficiently large 
detuning these can normally be neglected. 

All these terms conserve individually the z-projection of the total angular 
momentum of light and atoms, e.g. the S-j\ term can change a cr+ photon 
into a er_ photon (changing the light angular momentum along z by —2h while 
the atoms receive 2h mediated by the atomic raising operator j 2 + . The total 
angular momentum must have its ^-projection invariant since the physical 
system is axially symmetric around the direction of light propagation (the 
z-axis) . 

For us the term proportional to a\ is useful and relevant. This term represents 
the QND interaction. The higher order terms proportional to a 2 create some 
minor problems which can be minimized with large detuning (this will be 
discussed further in sections IB. II and IB.2|) . For a detailed treatment of th e 
higher order - atomic alignment - effects we refer to ( Kuprivanov et all 2005h . 



The zeroth order term proportional to ao produces an overall phase shift and 
can be omitted. 



2.4 PROPAGATION EQUATIONS 



The Hamiltonian (JJJ) is a very convenient starting point for many calcula- 
tions and we now show the procedure to derive equations of motion. For the 



spin operators we need the Heisenberg evolution dji/dt 



ih 



Ji 



H 



and for 



the Stokes operators the Maxwell-Bloch equations {d/dt + c- d/dz)Si(z,t) 



ih 



Si(z,t), H mi , see ([JulsgaardL 120031 ) . For the latter we will neglect retar- 
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dation effects, i.e. we do not calculate dynamics on the time scale L/c of 
propagation across the sample. This is equivalent to setting the speed of light 
to infinity and hence leaving out the d/dt term. If we consider only the term 
proportional to a\ and neglect the other we find 



dj x {z,t) cry A 2 d 

l2 



and 



dj z (z,t) 



dt 



0. 



d - 70 A 2 

— SJz, t) = H — aiSJz, t)j z (z, t), 

dz 8A2vr yV ; v y 

J^OM) = -^^^(M^OM), (10) 



We observe that j z (z, t) and S z (z,t) are individually conserved during the 
interaction which is also apparent from the Hamiltonian (J2J) since these op- 
erators commute with the ai-term. The effect of the interaction is that the 
spin will rotate around the z-axis with an amount proportional to S z , and the 
Stokes vector will rotate around the z-axis proportionally with %. 

Below we assume that these rotations are small and that the dominant clas- 
sical (mean) polarization vector of light and the orientation vector of the 
collective atomic spin stay oriented along the x-direction after the interaction. 
This turns out to be a very good approximation for experimentally attainable 
values of the interaction strength. Under this assumption the first line of the 
systems Q and (fTUjl. respectively, can be omitted. Furthermore, given the 
QND structure of the remaining equations, we can sum over the individual 
atomic spins and obtain the same equation for the collective spin variables ((TJ . 
In our continuous notation we have = J L %{z,t)pAdz. As for the light 
operators we concentrate on the in- and out-going parts only. Hence we define 
Sf 1 = cSi(z = 0,t) and S° nt = cSi(z = L,t). The multiplication by the speed 
of light c turns the normalization into photons per unit time instead of per 
unit length. With the assumption of small rotation angles, integrating Eqs. Q 
and (|r0|) over space from z = to z = L leads to the following very important 
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equations: 



Of 



S™(t)+aS x J z {t), 
aJ x S?(t), 



J 2 (t) =0, 



(11) 
(12) 

(13) 
(14) 



where a 



A 2 



- iA2^ ai - l n anc ^ ou ^ re f er to light before and after passing the 
atomic sample, respectively. 



We note from Eqs. (|11|) and (J 14)) that in the case of a large interaction strength 
(i.e. if aS x J z dominates S^ 1 ) a measurement on S° ut amounts to a QND mea- 
surement of J z . Using o ff- resonant light for QND measurements of sp ins has 
also been discussed in ( Kuzmich et al. . 19981 : Takahashi et al. . 1999f ). Equa- 
tion ()13j) implies that a part of the state of light is also mapped onto the atoms 
- we denote this as back action. This opens up the possibility of using this sort 
of system for quantum memory which will be the topic of sections 13 . 21 and f5 . 41 



2.5 THE ROTATING FRAME 



In the experiment a constant and homogeneous magnetic field is added in the 
x-direction. We discuss the experimental reason for this below. For our model- 
ing, the magnetic field adds a term Hb = £IJ X to the Hamiltonian. This makes 
the transverse spin components precess at the Larmor frequency fi depending 
on the strength of the field. Introducing the rotating frame coordinates: 



V 

v 

j: 



cos(fit) sin (fit) 
- sin (fit) cos (fit) 




(15) 



we can easily show that Eqs. (JTT|) - (|14j) transform into: 

ST'it) = S*(t) + aS x [j' y {t) sin(fit) + j' z (t) cos(fit 

aJ x S™(t) cos(fit), 



y 

^<out 

J^(t) = aJ x Sf{t) sin(fit). 



(16) 
(17) 

(18) 
(19) 



Thus, the atomic imprint on the light is encoded in the fi-sideband instead 
of at the carrier frequency. The primary motivation for adding the magnetic 
field is the fact that lasers are generally a lot more quiet at high sideband 



11 



frequencies compared to the carrier. A measurement without a magnetic field 
will be a DC measurement and the technical noise would dominate the subtle 
quantum signal. Also, as the measurement time is longer than 1/Q Eq. 
enables us to access both J' and J' at the same time. We are of course not 
allowed to perform non-destructive measurements on these two operators si- 
multaneously since they are non-commuting. This is also reflected by the fact 
that neither J y nor J z are constant in Eqs. (fTHj) and (fTT?|) . Below we shall con- 
sider two atomic samples and we will see that a QND-type interaction can be 
regained in this setting. 



2.6 TWO OPPOSITELY ORIENTED SPIN SAMPLES 



Inspired by the above we now assume that we have two atomic samples with 
oriented spins such that J x \ = —J x2 = J x - We re-express the equations of 
motion (|T^|) -()19 p for two samples in a way which is much more convenient for 
the understanding of our entanglement creation and verification procedure. 

For two atomic samples we write equations of motion: 

S*(t) + aS x ([J' yl (t) + J' y2 {t)\ Bin(nt) (20) 
+[J' zl (t)+J' z2 (t)]cos(Qt)), 

a(J x i + J x2 )S^(t)cos(nt) = 0, (21) 
a(J x i + J x2 )S?(t)sm(nt) =0. (22) 

The fact that the sums J' yl {t) + J' y2 {t) and J' zl {t) + J' z2 (t) have zero time 
derivative relies on the assumption of opposite spins of equal magnitude. The 
constancy of these terms together with Eq. (j2*Uj) allows us to perform QND 
measurements on the two sums. We note that each of the sums can be accessed 
by considering the two operators 

J o S° ut oos(nt)cfc = ^ S™cos(nt)dt+^(J' xl (t)+J' z2 (t)), (23) 

T T r T^ 

f Q 5° ut mn(Qt)dt = I S y «sm(nt)dt + ^(J yl (t) + J y2 (t)). (24) 

We have used the fact that J T cos 2 (ttt)dt w J T sm 2 (£lt)dt m T/2 and that 
Jq cos(fit) sin(Qt)dt ~ 0. Each of the operators on the left hand side can be 
measured simultaneously by making a ^-measurement and multiplying the 
photo-current by cos(fit) or sin(fit) followed by integration over the duration 
T. The possibility to gain information about J' y i(t) + J' y2 {t) and J' zl {t) + J' z2 (t) 
enables us to generate entangled states, the topic of sections 13.11 and 15.31 At 
the same time we must lose information about some other physical variable. 
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This is indeed true, the conjugate variables to these sums are j' z2 (t) — J' z i(t) 
and Jyi(t) — J'yzit)-, respectively. These have the time evolution 



§- t (J' y i(t) - .%,(*)) = 2aJ x S?(t) cos(fit), (25) 
§- t 0'At) - J'At)) = 2aJ x S?(t) frin(ni). (26) 



We see how noise from the input S^-variable is piling up in the difference 
components while we are allowed to learn about the sum components via S y 
measurements. The above equations clearly describe the physical ingredients 
for light and atoms. However, when transfer of a quantum state between two 
very different systems - light and atoms - is concerned, an isomorphic formu- 
lation, such as rendered by canonical variables, is desirable. In Sec. 0] we will 
re- write the theory in this more convenient language. 



3 Quantum Information Protocols 



In this chapter we introduce in detail two main protocols of this paper: de- 
terministic generation of entanglement of two macroscopic objects and deter- 
ministic quantum state mapping from one system onto another - the so-called 
direct mapping protocol. The latter is then applied as the quantum memory 
for light protocol. 

In order to simplify the description of the protocols we re-write equations of 
motion of the previous section in the language of canonical variables. This also 
allows for unified treatment of a single ensemble without- and two ensembles, 
with a constant magnetic field. For a single sample we define: 



^a - Jy Pa - Jz 

^ As ~ vX> As ~ v^*' ^7) 



For two samples we get two sets of canonical variables by defining the atomic 
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operators: 



*A. = J j^, (28a) 

p = ^1+^2 (2gb) 

*"=:w- (28c) 

P A2 = J -^k. (28d) 



2 /" T - 

' S y (t)cos(nt)dt, (29a) 
S z (t) cos(nt)dt, (29b) 
sin(fit)tZt, (29c) 
sin(Ot)d£. (29d) 
Each pair of X, P operators satisfy the usual commutation relation, e.g. we 




have 



i. Equations f)llti!4|) and (j20H22j) now translate into 

X°f = X™ + «Pg, (30a) 

Pu = Pu, ( 30b ) 
XXf = XZ + kPH (30c) 

Pm = PZ (30d) 

where we recall i = 1,2, s refer to the definitions above and not the two 
samples. Note, that in the case of two samples the two sets of interacting light 
and atomic operators are decoupled. The parameter describing the strength 
of light /matter- interactions is given by k = a\J 'J X S X T. The limit to strong 
coupling is around k ~ 1. This set of equations represents the starting point 
for numerous applications, which in the context of this paper are implemented 
in the double sample setup. 



3.1 ENTANGLEMENT - TWO MODE SQUEEZING PROTOCOL 



It follows from Eqs. ()30a|) and ()30c|) that a sufficiently precise measurement of 
Xl of a light pulse transmitted through two atomic ensembles renders knowl- 
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edge about P A . The measurement will proje ct Pa into a two-mode squeezed 
or Einstein- Podolsky- Rosen entangled state ((Julsgaard et all l2004aT ). 



The neces sary and sufficient condition for such an entangled state has been 
derived bv lDuan et aD (|2000h . Demonstration of entanglement is thus reduced 
to the fulfillment of this criterion. The criterion, which is introduced shortly, 
utilizes the vacuum state noise of a canonical variable, which for atomic ensem- 
bles reduces to the coherent spin state fluctuations. Hence, these fluctuations, 
which are also called projection noise, form an extremely important bench- 
mark in experiments with atomic ensembles. 



3.1.1 The Coherent Spin State 

To create entangled or squeezed states one has to generate states which exhibit 
less fluctuations than all equivalent classical states. The boundary occurs at 
the coherent spin state (CSS) in which all spins are independent realizations 
of a single spin oriented along a specific direction. The characteristics of this 
state are discussed in more detail in Sec. 15.11 so for now we need only to 
quantify the role of the coherent spin state as a boundary between classical 
and purely quantum mechanical states. From the canonical atomic operators 
a Heisenberg uncertainty relation can be formed: 



Var(X Ai )Var(P Ai ) > 



[X Ai ,P 



Ai 



1 

4' 



(31) 



For the coherent spin state the two variances are both equal to one half, 
thus confirming that it is the classical state with the least possible noise. 
For a state to be squeezed, it has to have less noise in one of the quadratures. 
Since the Heisenberg uncertainty relation still has to be fulfilled it follows that 
the other quadrature has to exhibit excessive fluctuations, it is anti-squeezed. 
Entanglement is the non-local interconnection of two systems, such that it is 
impossible to write the total density matrix as a product of density matrices 
for each system. It has been shown that for effectively continuous variable 
systems such as the ones described s o far in this paper the necessary and 



sufficient criterion for entanglement is ( Duan et al. . 2000) 



Var(P A1 ) + Var(P A2 ) < 1. 



(32) 



3.1.2 Entanglement Generation and Verification 

We now turn to the actual understanding of entanglement generation and ver- 
ification. We concentrate here on generation of an entangled state conditioned 
on the result of a measurement. Generation of an unconditional entangled 
state with the help of feedback is described in the experimental section below. 
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For generation of conditional entanglement we perform the following steps: 
First the atoms are prepared in the oppositely oriented coherent states cor- 
responding to creating the vacuum states of the two modes (-Xai,Pai) and 
(Xa2,-Pa2)- Next a pulse of light called the entangling pulse is sent through 
the atoms and we measure the two operators Xf^ and X£!f with outcomes A\ 
and Bi, respectively. These results bear information about the atomic opera- 
tors Pai and Pa2 and hence we reduce variances Var(PAi) and Var(pA2). To 
prove we have an entangled state we must confirm that the variances of Pai 
and Pa2 fulfill the criterion (J32j) . That is we need to know the mean values of 
Pai and Pa2 with a total precision better than unity. To demonstrate that, 
we send a second verifying pulse through the atomic samples again measuring 
X^ and X£!f with outcomes A 2 and B2. Now it is a matter of comparing A\ 
with A 2 and Pi with P 2 . If the results are sufficiently close the state created 
by the first pulse was entangled. 

Let us be more quantitative. The interaction (|3Ua|) mapping the atomic oper- 
ators Paj to field operators X Li is very useful for large k and useless if k -C 1- 
We will describe in detail the role of k for all values. To this end we first 
describe the natural way to determine k experimentally. If we repeatedly per- 
form the first two steps of the measurement cycle, i.e. prepare coherent states 
of the atomic spins, send in the first measurement pulse, and record outcomes 
Ai and Pi, we may deduce the statistical properties of the measurements. 
Theoretically we expect from ()30a|) 

1 K 2 

(A 1 ) = (Pi) = and Var(Ai) = Var(Pi) = - + — . (33) 

The first term in the variances is the shot noise (SN) of light. This can be 
measured in the absence of the interaction where k = 0. The quantum nature 
of the shot noise level is confirmed by checking the linear scaling with photon 
number of the pulse. The second term arises from the projection noise (PN) 
of atoms. Hence, we may calibrate k 2 to be the ratio k 2 = PN/SN of atomic 
projection noise to shot noise of light. Theoretically k 2 has linear scaling k 2 = 
aJ x S x T with a macroscopic spin J x that must be confirmed in the experiment 
(see Sec. 15. 

Next we describe how to deduce the statistical properties of the state created 
by the entangling pulse. Based on the measurement results A\ and Pi of this 
pulse we must predict the mean value of the second measurement outcome. If 
k — ► 00 we ought to trust the first measurement completely since the initial 
noise of Xf^ is negligible, i.e. (A 2 ) = A x and (P 2 ) = Pi. On the other hand, 
if k = we know that atoms must still be in the vacuum state such that 
(A 2 ) = (B 2 ) = 0. It is natural to take in general (A 2 ) = aA 1 and (P 2 ) = a.B\. 
We need not know a theoretical value for a. The actual value can be deduced 
from the data. If we repeat the measurement cycle N times with outcomes 
Af, Pf } , Af, and Bf, the correct a is found by minimizing the conditional 
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variance 

Vax{A 2 \A x ) + Vax{B 2 \B x ) = 

TiV^T^l^ aA i ) +(S^ ; -a J B| 7 ) J. 

i 

In order to deduce whether we fulfill the entanglement criterion ()32|) we com- 
pare the above to our expectation from (j30ap . For the verifying pulse we get 



(*S? - (*u )) ) = ( ( 



' -iMn,2nd i 



pent 



P 



cnt \ 



Ai / 



- + ^ 2 Var(Plf ) 



(35) 



where X L °' 2nd refers to the incoming light of the verifying pulse which has zero 
mean. PJg* refers to the atoms after being entangled. We see that the practical 
entanglement criterion becomes 



Var(A 2 |Ai) + Var(P 2 |Pi 



1 



K 



(Var(PXf ) + Var(P^)) 



< 1 + k 2 = Var(Ai) + Var(Pi). 



(36) 



Simply stated, we must predict the outcomes A 2 and B 2 with a precision 
better than the statistical spreading of the outcomes A x and B x with the 
additional constraint that A x and Pi are outcomes of quantum noise limited 
measurements. 



3.1.3 Theoretical Entanglement Modeling 

We described above the experimental procedure for generating and verifying 
the entangled states. Here we present a simple way to derive what we expect 
for the mean values (i.e. the a-parameter) and for the variances Var(P^). 

We calculate directly the expected conditional variance of A 2 based on A x : 

-{>out,2nd -{>out,l 
X L1 ~ aX Ll 



(x L f nd - aXft™ + k [PS - aPlf]) 2 ) (37) 



=^(l + « 2 + K 2 (l-af). 

In the second step we assume that a perfect QND measurement without any 
decoherence is performed, i.e. P^f = P^. By taking the derivative with re- 
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spect to a we obtain the theoretical minimum 



Var(A 2 |A x ) + Var(5 2 |fi 1 ) = 1 + 1 

1 + «r 



Var(PXi t ) + Var(P^) 
obtained with the a-parameter 



l + /t 2 



a = j^. (39) 

It is interesting that, in principle, any value of k will lead to creation of 
entanglement. The reason for this is our prior knowledge to the entangling 
pulse. Here the atoms are in a coherent state which is as well defined in terms 
of variances as possible for separable states. We only need an "infinitesimal" 
extra knowledge about the spin state to go into the entangled regime. 

It is also interesting to see what happens to the conjugate variables Xm in 
the entangling process. This is governed by Eq. ()30c|) . We do not perform 
measurements of the light operator so all we know is that both and 
Pl" are in the vacuum state. Hence Var(X^f) = (1 + k 2 )/2 and we preserve 
the minimum uncertainty relation Var(X| 1 J t )Var(P^° t ) = 1/4. 



3.1.4 Entanglement Model With Decoherence 

Practically our spin states decohere between the light pulses and also in the 
presence of the light. We model this decoherence naively by attributing the 
entire effect to the time interval between the two pulses, i.e. we assume there 
is no decoherence in presence of the light but a larger decoherence between 
the pulses. We may then perform an analysis in complete analogy with the 
above with the only difference that P^f = PPai + \/l ~ P 2 Vp where V p is a 
vacuum operator admixed such that (3 = corresponds to a complete decay 
to the vacuum state and [3=1 corresponds to no decoherence. Completing 
the analysis we find the theoretical conditional variances 



1 _L (1 _ fp\ 2 

Wai(A 2 \A 1 ) + VarOBalBx) = 1 + ' 



=► Var(P^) + Var(PX™ 
obtained with a-parameter 



1 + K 

■ ut \ - ,., n< ni > 1 + (1 — )k 



2 



(40) 



« = T^Ti- (41) 



[3k 2 

l + K 1 

In the limit /3 — > 1 these results agree with (|38|) and (|3*^j) . For (3 — > we have 
a — > (outcomes A\ and E>i are useless) and the variance approaches that of 
the vacuum state which is a separable state. 
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3.1.5 Gaussian State Modeling 



Following the extensive study of the evolution of Gauss i an states during arbi- 
trary interactions and measurements ( Giedke and Cirac , 20021: Eisert and PlenioL 
20031) the development of spin squeezing in a single atomic s ample and entan- 



glem e nt in two samples were treate d in detail theoretically in ([Madsen and Molmer , 
l2004t ISherson and Molmer . 2005). Arbitrary Gaussian states of n canonical 
modes are fully characterized by a 2n x 1 vector, v, describing the mean val- 
ues and a 2n x 2n matrix, 7, describing the correlations within the atomic and 
light systems and the cross-correlations between these. In this formalism an 
interaction between light and atoms is governed by: v — > Sv and 7 — > S'jS 7 '. 
For the Faraday interaction, the coefficients of the matrix S are easily deter- 
mined from Eqs. (|30j). Decoherence is easily included but most importantly, 
there is an explicit expression for the state of the remaining modes after an 
arbitrary homodyne measurement on a number of the modes is performed. If 
the incoming pulse of light is split into a large number of segments and the in- 
teraction with atoms and subsequent measurement of each segment is treated 
sequentially differential equations for v(t) and j(t) can be obtained in the 
limit of infinitesimal segment durations. The former is a stochastic differential 
equation determined by the outcome of the measurement, whereas the latter 
is d eterministic, although non- l inear because of the measur ement dynamics. 
In ( Madsen and Molmerl . l2004t ISherson and Molmer . 2005| ) such differential 
equations are solved giving the time resolved dynamics in the presence of e.g. 
atomic decoherence and light losses. In this way, e.g. analytic expressions for 
the optimal degree of spin squeezing and the degree of entanglement at arbi- 
trary rotation frequencies (the verification and generalization of Eq. (|38j) ) can 
be obtained easily. 



3.2 QUANTUM MEMORY 



For complete quantum memory we require 1) that the light state to be stored 
is supplied by a third party in an unknown state, 2) that this state is mapped 
onto an atomic state with a fidelity higher than the best classical fidelity, and 
finally 3) that the stored state can be retrieved from memory. As is described in 



Sec. I5.4l the first two criteria have been met experimentally in (lJulsgaard et al 



2004af l whereas the last one still remains an unsolved experimental challenge. 



A recently developed experimentally feasible protocol for retrieval now exists 
and is discussed in Sec. 13.2.31 
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3.2.1 Direct Mapping Protocol 



In Eqs. (l30aH30c)) one of the light variables is mapped onto one of the atomic 
variables. This represents a natural starting point for a quantum memory 
protocol in which the entire light mode described by the two non-commuting 
variables X m and P m is faithfully s tored. In the so-called "direct mapping 



protocol" of (jJulsgaard et all l2004a[ ) the mapping is completed by measuring 
the remaining light quadrature x out = x m + kP" 1 and feeding the result back 
into the atomic Xa with a gain of g: 

X° A ut = XX 1 + nPi n , (42a) 

pout , = pout _ ^out = pm (1 _ Kg) _ ^in ^ 

If k — g — 1 and the initial atomic state is assumed to be a coherent state 
with zero mean value the mean values of both light variables will be stored 
faithfully in the atoms. Although the initial atomic state has zero mean it 
is a quantum mechanically fluctuating state and any uncanceled atomic part 
will increase the noise of the final state and thus degrade the mapping perfor- 
mance. Although this protocol works for any state, in the following we discuss 
storage of coherent states of light, i.e. vacuum states which are displaced by 
an unknown amount in phase space. For the storage of an arbitrary coherent 
light state the remaining X A contribution limits the storage fidelity to 82%. 
This can be remedied by initially squeezing the atomic state, in which case 
100 % fidelity can be reached in the limit of infinite squeezing. 



3.2.2 Mapping With Decoherence 

Just as in the case of entanglement generation the spin states decohere. Again 
we can model this by a beam splitter type admixture of vacuum components 
right after the passage of the first light pulse. We can furthermore model light 
damping (e.g. reflection losses) in a similar way to obtain: 

XX ut - + kP£) + y/l-PVxA, (43a) 

P£* - (0 - gKy/c)P% - gy/tX* + sfl - [3 2 V PA - s/hVxl • (43b) 

We see that and X™ are mapped with gains g' BA = (3k and g' F = g^fC, 
respectively. The variances can be calculated easily to be: 

Var(X out ) = l + sg A , (44a) 
Var(P out ) = l + f + - 2g' F g' BA . (44b) 
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3.2.3 Quantum Memory Retrieval 



As mentioned in ( Julseaard et al. . 2004eJ ) the stored state can in principle be 



retrieved by inverting the roles of light and atoms in the direct mapping pro- 
tocol. This would involve first an interaction between a read-out light beam 
with the atomic sample acting as a storage medium. According to Eq. ()30a|) 
this would map Pa onto the light. Next has to be measured and feed- 
back applied to the read-out beam according to the result of the measure- 
ment. However, since the atomic measurement requires a certain time during 
which the read-out pulse propagates at the speed of light, the feedback is 
only prac tical for pulse durations shorter than a microsecond. In the exper- 



iments of lJulsgaard et al.l (|2004ar ) pulses of millisecond duration (~ 300km) 
are required in order to obtain a sufficiently high interaction strength, and 
the inverse direct mapping protocol is thus infeasible for this experimental 
realization. 

Several years ago a retrieval scheme, which did not involve measurements, 
but instead used two orthog onal passages of the read-out pulse was proposed 



(jKuzmich and Polzikl . 120031 ) . In each passage one of the atomic quadratures is 
mapped onto the light pulse, and in this way retrieval fidelities of up to 82% 
can be achieved with coherent readout light and up to 100% with perfectly 
squeezed readout light. Unfortunately in order to preserve the QND nature 
of each of the two interactions the light has to pass entirely through the 
atomic medium before proceeding to the second passage. This again renders 
the protocol inapplicable to all setups requiring "long" pulses. 

Recently this problem was eliminated by solving the complex dynamics arising 
when a light beam passes through the atomic m edium along two orthogonal 
directions simultaneously (|Sherson et al. . 2005b| ). This protocol works with 



arbitrary pulse durations and the fidelity of this two-pass protocol has been 
calculated both for a coherent input state and for a light qubit. It was also 
shown that if the light is reflected back after the second passage, thus com- 
pleting four passages, and a time dependent interaction strength is applied 
perfect retrieval ca n be achieved withou t requiring squeezed initial states. In 
a related proposal ( Fiurasek et al. . 20051 ) retrieval is achieved by sending two 



different beams through the atomic samples simultaneously, each in a direction 
orthogonal to the other. The fidelity of retrieval for coherent states including 
realistic light losses is calculated and shown to exceed the classical bound. 



3.2.4 Alternative Quantum Memory Proposals 

In (|Kuzmich and Polzikl . I'M)! various protocols which can be implemented 
using the QND-Faraday interaction are described in detail, including atom- 
atom teleportation. Instead of reviewing these, we would like to mention 
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some recent theoretical proposals for improving the quantum memory per- 
formance. All of these break with the simple interactio n scheme either by in- 
troducing a nqn-QND interaction (as was also done in (|Sherson et all l2005bt 



Fiurasek et al. . 2005| )) or by exploiting the fact that the atoms involved are 



not simple spin- 1/2 atoms and thus contain more than two magnetic sublevels. 



In f Opatrnv and Fiurasek . 2005| ) it is proposed to exploit additional atomic 



coherences to enhance the capacity of an atomic quantum memory. Remember 
that during the usual QND Faraday interaction Am = 1 atomic coherences 
are coupled to the sideband of the light at the Larmor frequency. If in addition 
circularly polarized classical fields are present during interaction coupling be- 
tween Am = 2 coherences and a light sideband at twice the Larmor frequency 
are also created. In this way an additional quantum channel is added to the 
quantum memory, thus enhancing its capacity. By appropriately tuning the 
detuning of the additional coupling field two-mode squeezer and beam split- 
ter Hamiltonians can also be realized. With the latter, a quantum memory is 
realized in a single passage without a requirement of measurements or prior 
squeezing of atoms or light. 



In (|ODatrnvi . EoOfih the author proposes to create a single cell atomic memory 



by optically pumping the sample into an incoherent mixture of the two extreme 
magnetic sublevels (m = —4 and m = 4 for the F = 4 hyperfine level of 
cesium). This is done with suitable linearly polarized light. When sending 
a light pulse through the atomic sample it will interact with both extreme 
coherences simultaneously, creating the usual QND-Faraday interaction. Apart 
from the prospect of making the memory more compact, this proposal would 
also substantially decrease the reflection losses associated with the many air- 
glass transitions in a two-cell setup. A main problem in this approach is that in 
the presence of the probe light the two coherences will experience different AC- 
Stark shifts and thus different phase evolutions. For relevant pulse durations 
this effect is in fact significant. The author proposes to solve this by either 
decreasing the bias magnetic field or by introducing an additional AC-Stark 
shift by adding an light field with appropriate polarization and detuning. 

An atomic memory could of course also be implemented using light-atom tele- 
port at ion. This could be achieved by first sending an auxiliary pulse through 
the two oppositely oriented atomic samples. This would entangle the light and 
atomic systems via the QND-Faraday interaction. After this, the auxiliary 
light beam should be mixed on a 50/50 beam splitter with the quantum light 
signal to be teleported. Measurements of S y and S z in the two output ports 
respectively and subsequent feedback to the atomic system would complete 
the teleportation. Unfortunately the achievable fidelity is li mited to 67% with 



the us ual QND-interaction. It was therefore proposed in ([Hammerer et al 



2005af ) to replace the two atomic samples in the protocol just sketched by a 



single sample, still with a constant bias magnetic field. This changes the in- 
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Figure 2. (a) A photographic view of a typical experimental setup. Atomic vapor 
cells are placed inside the cylindrical magnetic shields. The pumping beams are 
indicated with dashed arrows and the path of the quantum probe field is marked 
with the solid arrows, (b) A schematic view of the setup. The probe pulses reach a 
detection system measuring S y (t). The photo current is sent to a lock- in amplifier 
which singles out the sin(f^) and cos(Ot) parts as we discuss in Sec. I4.2I These 
are integrated and stored in a PC. Typically, the pulse sequence consists of (1) a 
pumping pulse, (2) and (3) two laser pulses for quantum manipulation and detection 
of the atomic states. A very small portion of the probe field is sent through each 
sample in the x-direction to measure the magnitude J x of the macroscopic spins by 
Faraday rotation measurements. 

teraction dramatically. Both atomic quadratures are transferred to the light, 
albeit at the cost of additional new modes of light. The resulting multi-mode 
light-atom entanglement enables teleportation with up to 77% fidelity. 



4 Experimental Methods 



In this section we describe the typical setup for our experiments. It is centered 
around two glass cells filled with Cs vapor at room temperature placed in 
two separate magnetically insulating shields with a bias magnetic field inside. 
Additional coils are used to apply an rf magnetic field with the frequency 
equal to the Larmor frequency of the bias field. Atoms in the cells are optically 
pumped. The setup, therefore, is similar to that of a classical magneto-optical 
resonance experiment. 

The typical experimental setup is shown in Fig. |2j A Verdi V8 pumped 
Thsapphire laser delivers what we call the probe field. This is used for quan- 
tum manipulation of the atomic samples and also for detecting the macroscopic 
spin J x . Diode lasers provide optical pumping fields for creating highly polar- 
ized spin states. The cw lasers are modulated by AOMs or EOMs to achieve 
the desired pulse shapes. Atoms are contained in vapor cells and placed in 
stable magnetic fields. The magnetic field homogeneity must be of order 10~ 3 
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Figure 3. The setup for measuring S y . A strong classical pulse is mixed on a polariz- 
ing beam splitter with the quantum field which can be in a vacuum state, as in the 
entanglement experiment. Measurement is performed in the 45°- and — 45°-basis 
by two detectors. Half of the difference of the two photo currents is S y . An op- 
tional A/4-retardation plate turns this measurement into the measurement of S z 
component. 



across the vapor cell volume. 



4.1 PARAFFIN COATED VAPOR CELLS 



In our experiments the atomic samples are contained in a paraffin coated vapor 
cell. The coating prevents depolarization of the spin state when atoms hit the 
walls. We have measured spin coherence times exceeding 40ms and spin po- 
larizations exceeding 99% by the methods described in Sec. 14.3.11 The atomic 
density and thereby the macroscopic spin J x can be controlled by heating or 
cooling the vapor cell. In order to achieve a stable vapor density, tempera- 
ture gradients across the cell should be avoided. Temperature control by air 
flow is a convenient solution. Metal heating/cooling elements causes severe 
problems since the atoms are disturbed by rando m magnetic field s created 
by thermal currents even if aluminum is used, see (l.lulsgaardl l2003h. Further 
information on paraffin coated cells can be found in (jBouchiat and Brossel 
196i lAlexandrov et all Fiji liffll. 
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4.2 DETECTION OF POLARIZATION STATES 



The Stokes parameters are measured with low noise photo detectors. We use 
high quantum efficiency photo diodes and home made amplifiers characterized 
by negligible electronic noise compared to the shot noise of light at optical 
power higher than lmW. In Fig. |3] we depict how the S^-component of light 
is measured. 

The differential photo current i(t) from the two detectors corresponds to a re- 
alization of the measurement of S y . By passing i(t) through a lock-in amplifier 
we can detect the sine and cosine components at the Larmor precession fre- 
quency Q. Practically, the current i(t) is multiplied by cos(f2t) and sin(f2t) and 
integrated over time for this purpose. According to (j2HH-d), with appropriate 
scaling, this exactly corresponds to measuring the Xli and Xl2 components 
of light. For ^-detection we would measure Pli and Pl2- 



4.3 MAGNETIC FIELDS 



As was discussed in Sec. 12.51 external magnetic fields are added to control the 
dynamics of the atomic samples. The interaction of atoms with a magnetic 
field is governed by the Hamiltonian 

# mag = ^ B J-B + 0(P 2 ) (45) 

We stress here that J is the total angular momentum of the atom including the 
nuclear spin. For the F = 4 ground state of cesium gp « 1/4. The second term 
0(B 2 ) reminds us that the above linear equation is only approximately true. 
When the magnetic energy becomes comparable to the hyperfine splitting of 
the ground state the response is non-linear. We comment on this below. 



Adding a constant bias magnetic field B x in the x-direction leads to the equa- 
tions of motion (JT5J) with Larmor precession at the frequency = 9F^BB x /h. 
If we furthermore add an RF magnetic field at frequency Q along the y- 
direction such that 



B C xt = B x e x + (B c cos(fit + 4>) + B s sin (fit + <p)) e y (46) 

with constants B c and B s we may derive for the rotating frame coordinates 
J' y and J' z of (USD that 

— ^— = — lo s sin(fi L £) sin(£7t + 4>)J X , ^ ^ - = cos(il^t) cos(fit + 4>)J X , 

(47) 
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with u CyS = gpii^Bc^jU. Choosing the phase and the frequency of the RF-drive 
such that = and fl = we obtain: 

j. (t ) = j.(o), 8 Jm = -^k, a|w = _^. (48) 

w w ' dt 2 dt 2 v ' 

These equations are valid for interaction times T such that u c T, uj s T <C 1 <C 
VtT . We see that with pulses of RF-magnetic fields we are able to change the 
spin components J' y and J' z by an amount controlled by the sine and cosine 
components B s and B c . This has several experimental applications, which are 
discussed below. 



4.3.1 Characterizing the Spin State with the Magneto-optical Resonance 
Method 

Equations (j20J) and (JI7J) describe the Magneto-Optical Resonance method 
(MORS) which we use extensively for the spin state charact erization. Appli 



cation o f MORS to our experiments are described in detail in iJulseaard et al. 
(2004iJ. Within this method the RF frequency is scanned around and the 
oscillating transverse spin components are probed via oscillating polarization 
rotation of the optical probe. In order to quantitatively explain the MORS 
signal as the RF- frequency is scanned across Q^, we need to turn back to the 
second order term mentioned in Eq. ()45|) . The transverse spin components j y 
and % can be expressed in terms of coherences <3" m ,m±i m the following way: 



h = 2 51 y/ F ( F + 1 ) - m ( m + l)(0Vn+l,m + CTm,m+l), 



(49) 



3* = 2i^l2 y F ( F+ X ) ~ m ( m + 1 )(^m+l,m - ^m,m+l)- 



In the absence of the second order term in ()45|) the energy separation hQ^ 
between states \m) and \m + 1) is the same for all m and all terms <3" m +i,m 
have the same resonant frequency. The second order term in ([45)1 . however, 
makes the frequency of the coherences oVn+i.m slightly different. It can be 
shown, that the frequency difference u>qz between a m ^ m +i and a m -x, m is wqz = 
2f2L/o;hfs where tUhfs = 27r ■ 9.1926GHz is the hyperfine splitting of the Cesium 
ground state. We typically have = 2tt ■ 322kHz and the effect is small but 
detectable. 

In the special case that the amplitude and frequency of the driving RF-field 
vary slowly compared to the spin coherence time, the off-diagonal coherences 
follow the diagonal populations adiabatically and we may write e.g. J y as: 



J y = Re 



^ [F(F + 1) — m(m + 1)] ■ e int 

const -— — — — j— [a m+ i >m+ i — (J m ,rn\ 



(50) 
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Figure 4. MORS signal, in arbitrary units (a.u.) of a poorly oriented sample. Each 
peak corresponds to a « + 1 Zeeman resonance. The height of each reso- 
nance is proportional to the population difference of the two relevant sublevels. 

where r m+l m are the FWHM linewidths giving an exponential e~ rt ^ 2 decay of 
each coherence. The Larmor frequency has been replaced by the individual 
coherence evolution frequencies fi m+l m . For J z we have to take the imaginary 
part. Two adjacent magnetic sublevels act as a two level system with the 
usual Lorentzian response to a driving RF field. Scanning the RF frequency 
we get eight Lorentzian peaks, the magnitudes of which will depend on the 
populations of the magnetic sublevels. The MORS signal is proportional to 
the square of the term in the square brackets in Eq. ()50|) . An example of such a 
signal can be seen in Fig.^J The application of this method is twofold. First, we 
gain information on the distribution of population among the different ground 
state magnetic sublevels. From this we infer that we are able to optically pump 
the atoms to such an extent that only the outermost coherence (m? = 4 <-» 
mp = 3) becomes significant. Second, measuring the widths of the resonances 
under different experimental conditions allows us to quantify the effect of 
different decoherence mechanisms as discussed in Sec. 15.21 

A second implication of the second order term is that the equations of motion 
djy/dt = —Qjz and dj z /dt = Qj y describing Larmor precession are too simple 
and must in principle be generalized. However, this is not necessary for our 
case with highly polarized samples. We have almost all the population in the 
mp = 4 and mp = —4 states for the two samples. Then there is effectively 
only one non- vanishing frequency component in the sums (J49|) . For the j x = 4 
sample the single term is 034 and for the j x = —4 sample the only term is 



4.3.2 Manipulating the Spin State 

Of course, as suggested by Eq. (f4~Sj). an appropriate choice of phase, strength 
and envelope function of the RF-field allows us to create an arbitrary mean 
value of either of the spins or any combination of these. This has two appli- 
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cations. First, it enables us to create a large classical mean value and observe 
the evolution of this state, thus constituting another mechanism for studying 
decoherence and calibrating the system as is discussed briefly at the end of 
Sec. IB.ll Second, we can actively feed back the result of a quantum probing 
of the spin state, thus creating a particular desired state. This is the keystone 
element for creation of deterministic entanglement, discussed in Sec. 15.3.2*1 and 
for the quantum mapping experiment discussed in Sec. 15.41 



5 Experimental Results 

5.1 PROJECTION NOISE LEVEL 

Since the Heisenberg uncertainty relation sets the starting point of all our cal- 
culations, one of the most important tasks in our experiments is the achieve- 
ment of quantum noise limited performance. Practically, it is also one of the 
most difficult tasks. When we detect polarization states of light we observe 
noise in the signals. After the light has passed the atomic samples, there is a 
contribution to this noise from the light itself and from the atomic spins. The 
noise contribution from atoms in the minimum uncertainty state (the coherent 
spin state) is called the projection noise. 

We discussed the ratio of the projection noise to the quantum noise of light 
(shot noise) already in Sec. 13.11 We found that theoretically this ratio should 
be 

k 2 = a 2 J x S x T. (51) 

The ratio k 2 came from the more general interaction equations (|3tM d). In the 
present section we discuss how to calibrate this projection noise level experi- 
mentally and how to predict the noise level from independent measurements. 

5.1.1 Measuring the Macroscopic Spin 

The ratio of projection noise to shot noise is proportional to the macroscopic 
spin J x . This linearity is the finger print of quantum noise and is essential to 
establish experimentally. To this end we need a good measure of J x . 

We measure J x by detecting polarization rotation of a linearly polarized laser 
field propagating through the atomic samples along the x-direction. To see 
what happens in this setting we consider Eqs. Q and (jlUj) . These equations 
assume propagation along the ^-direction so we assume the spin to be polarized 
along the ^-direction in the following. For linearly polarized light we have 
S z j = and the effect on the transverse spin components j Xl j y is negligible. 
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It can also be shown easily that the ao and ci2 terms of the Hamiltonian (|7|) 
play no role in this calculation. We are left with the evolution of S x and S y 
according to (|l(Jj) and after integration over the sample we find 



ST" = SJT cos(2ft,) - 5; n sin(2^ F ), 

S° y ut = S x n sm{29 F ) + S y n cos{29 F ), 1 j 

where "in" refers to the polarization state before the sample at z = and 
"out" refers to the state after the sample at z = L. The angle 9-p is given by 
(in radians) 

If a linearly polarized beam of light is rotated by the angle 9, the Stokes vector 
is rotated by 29. Thus, in the above, 9 F is the polarization rotation caused by 
the spin orientation along the direction of light propagation. We note that 
the angle 9 F depends on the density p of atoms times the length L that the 
light traverses. We wish to re-express this in terms of the macroscopic spin size 
Jz = N at {%) of the entire sample (remember we have the spins polarized along 
z in this discussion). To this end we observe that iV a t = pV = pA ce \\L where V 
is the vapor cell volume and A ce n is the area of the vapor cell transverse to the 
beam direction. This will conveniently be an effective area for cells that are 
not exactly box like. Turning back to the usual convention of spin polarization 
along the x-axis we then rewrite (J53|) as 



327Tv4 cell A 



5.1.2 Predicting the Projection Noise Level 

Now, let us return to the predicted ratio of projection to shot noise JSTJ. This 
prediction relies on Eqs. (|lllll4j) which are derived under the assumption that 
all atoms in the sample are seen by the laser beam which has a cross sectional 
area A. In experiments the laser beam does not intersect all the atoms. In 
App. El we show that the random motion of atoms in and out of the beam 
modifies the expected variance of the transverse spin components J y and J z by 
statistical effects from the usual J x /2 to p 2 (l + a 2 )J x /2 where p = A/A cea is 
the mean time of an atom inside the laser beam and a 2 is the relative variance 
of p. We furthermore present a simple model for p and a and show that the 
atomic motion will act as an effective source of decoherence between two probe 
pulses. We incorporate atomic motion into Eq. (joTj) by replacing A with A ce u 



29 



2 4 6 8 10 12 14 16 

Faraday angle [deg] 

Figure 5. Measured atomic noise rel- 
ative to shot noise of light. The lin- 
earity is a clear signature of the pro- 
jection noise limitation. The slope 
K exp = 0.104(2) • f?p should be com- 
pared to the theoretical value of 
K 2 h = 0.140 -9 F from Eq. (JHHJ) with 
a 2 = 0. 
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Figure 6. Slope of measured k? 
vs. 6 normalized to the experimen- 
tally predicted level (without the fac- 
tor 1 + a 2 ) vs. Tp ro be- The fit gives 
/4 p //c 2 h = l + 0.47(13)/T[ms]. 



in the factor a and multiplying the whole expression by 1 + a 2 . We then find 



K 2 = a 2 J x S x T-p 2 (l + a 2 ) = (— ^— ^ ai ) J X S X T(1 + a 2 ) 



8A cell A 2tt 



l + a 2 )7A 3 aiP-r-fj F 
327r 2 A ccn A^c ' 



(55) 



In the last step we replaced S x = <p/2 = P/2tkv = PX/AkTic where P is 
the optical power. We also inserted Eq. (p>4"j) to express k 2 as a function of # F . 
However, we must remember that the area A ce n in (p>4"j) refers to the transverse 
area for a beam propagating in the x-direction while the A ce n from the relation 
p = A/A cc n refers to the transverse area for a beam along z. Hence, the last 
step above is valid for a vapor cell of cubic symmetry only, but it can still be an 
irregularly shaped cell. In other cases the generalization is straightforward. We 
have reached an expression for k 2 in terms of convenient parameters from an 
experimental point of view. With 7/27r = 5.21MHz and A = 852. 3nm together 
with he we reach our final theoretical estimate for the projection to shot noise 
ratio expressed in convenient units: 



2 _ 56.4 ■ P[mW] ■ T[ms] ■ fl F [deg] ■ fll (A) ■ (1 + a 2 ) 

K th - a r 21 A rA/ru„l ' (5b) 



A cell [cm 2 ]-A[MHz] 



where cq(A) was denned in Eq. (JBJ). 
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5.1.3 Experimental Investigation 



Turning to experiment, in Fig. we see an example where the measured noise 
relative to the shot noise of light is plotted. The data are clearly linear. With 
A/2tt = 700MHz, T = 2.0ms, P = 4.5mW, and a 2 = for the moment, we 
predict a linear slope of 0.140 which is somewhat higher than the measured 
value of 0.104. Including the a 2 from atomic motion makes the discrepancy 
slightly worse but given the simplicity of the theoretical model we consider 
the agreement satisfactory. 

To test the scaling properties predicted in the atomic motion calculations, 
we fix the power P, detuning A and macroscopic spin size J but vary the 
probe duration T. The measured noise is plotted in Fig. relative to the 
prediction (JoTlj) with a 2 = 0. We see that as T is increased we do see a lower and 
lower noise level which corresponds to decreasing a 2 . The solid line in the figure 
represents a fit where a 2 = (0.47±0.13)/T[ms]. To compare this to the simple 
model described in (jA.4|) we estimate our beam diameter to be 1.6cm which 
gives A « 2.0cm; moreover, L = 3.0cm, vq = 13.7cm/ms (cesium at room 
temperature). For T = 1ms we get the prediction a 2 = 0.44. This is in very 
good agreement with the measured data, but this agreement must be viewed 
as fortuitous. As mentioned before, numerical simulations of atomic motion 
have shown that the variance estimate (|A.4j) is almost four times too high. 
The high experimental value must be attributed to the additional Doppler 
broadening effect. We also note the relatively high uncertainty of 0.13. But all 
together we have a qualitative understanding of the physics and a quantitative 
agreement within a few tens of percent. 



5.1.4 Thermal Spin Noise 

Another issue concerning the projection noise level is the question of thermal 
spin noise. For the establishment of the correct noise level we must be in the 
CSS with high precision. For the CSS the spin is completely polarized along 
the x-direction and Var(j y ) = Var(j^) = F/2 = 2 for the F = 4 ground state. 
As a very different example we may consider a completely unpolarized sample. 
We then have by symmetry Var^) = V&r(j y ) = Var(j z ) = (j 2 + j 2 + j 2 )/3 = 
F(F + l)/3 = 20/3. This is a factor of 10/3 higher and, even for fairly good 
polarization, the thermal noise may be significant. In our experiments with 
quantum information protocols we exceed a spin polarization of 99% which 
means that the thermal noise must be very small compared to the true projec- 
tion noise. The degree of spin po larization has been mea sured independently 
with methods similar to those in lJulsgaard et al. (|2004bh . 



A nice illustration of the fact that we get lower noise for the CSS than in the 
unpolarized case is given in Fig. [7| Experimentally, we perform measurements 
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y=0.342(1 7)(exp(-x/7.8(1 5))+2.05(9)(1 -exp(-x/7.8(1 5)))) 



Time [ms] 

Figure 7. Coherent state noise compared to the completely unpolarized spin noise. 
The data is taken with a vapor cell in which the spin life time is very short. The 
noise level increases on a time scale of roughly 8ms to the thermal equilibrium level. 
The increase in noise is consistent with predictions for the coherent and unpolarized 
spin states. 

on very poor vapor cells where the macroscopic spin life time is small. We op- 
tically pump the sample and wait for some variable delay time before probing 
the spin noise. For long times the spins will reach thermal equilibrium, where 
the noise of each atom in F = 4 contributes 20/3. The fraction of atoms in 
F = 4 is 9/16, the remaining 7/16 are in the F = 3 state and do not contribute 
because of the large detuning. Initially, all atoms are in F = 4 in the CSS and 
they each contribute the value 2 to the noise. Hence the measured noise must 
be on the form 

20 9 

Measured noise oc 2 • exp(— Ft) H (1 — exp(— Tt)). (57) 

3 16 

The predicted ratio of final to initial noise is thus 15/8 ~ 1.88. Experimen- 
tally we find the ratio 2.05 ± 0.09 which is consistent. To sum up, there is 
strong evidence that we really do create the CSS with the correct minimum 
uncertainty noise. 



5.1.5 Concluding Remarks on the Projection Noise Level 

Let us sum up the discussion of the projection noise level. To reach the quan- 
tum noise limited performance one should first of all observe the atomic noise 
grow linearly with the macroscopic spin size J x . An experimental example of 
this was shown in Fig. El The linearity of the noise basically arises from the fact 
that different atoms yield independent measurements when their spin state is 
detected. Technical noise sources from e.g. external electromagnetic fields cou- 
ple to all atoms and the effect on the noise variance would be quadratic. 
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However, linearity alone is not enough. An ensemble of independent and unpo- 
larized atoms would also show a linear increase in the spin noise when increas- 
ing the number of atoms. Since unpolarized atoms have larger noise variance 
than the 100% polarized atomic sample, we must know independently that 
the spin orientation is high. In our experiments the spin samples are polarized 
better than 99%. 

One may argue that the small fraction of atoms that are not in the completely 
polarized state could, in principle, form exotic-multi particle states with a very 
high variance of the detected spin noise. The results discussed in Sec. 15.1.41 
prove that this is not the case. 

Finally, as derived in App. |X] the atomic motion leads to an increased ratio 
of atomic to shot noise. Generally a large ratio of atomic projection noise to 
shot noise is good for the quantum information protocols. However, we do not 
gain anything by the increase of the atomic noise caused by atomic motion 
since there will be an accompanying increase in the decoherence rate. 

5.2 DECOHERENCE 

As mentioned earlier, all atoms are optically pumped into an extreme Zeeman 
sublevel with the quantization axis. A conventional way of categoriz- 

ing sources of decoherence is according to whether they affect the magnitude of 
the spin along this axis or merely along transverse directions. The appropriate 
life times of these are called Ti and T 2 and defined as: 

J x (t) = e- t/Tl J x (0) and J trans (t) = J tr ans(0) (58) 

in the absence of additional interactions. In the absence of the probe light 
the dominant processes are collisions with the walls and other atoms and, 
typically T\ w 300ms. Even though the probe detunings are quite high in our 
experiments (700-1200MHz) making the desired refractive Faraday interaction 
dominant by far, the small probability of absorption still reduces T\ by about 
a factor of 2 depending on probe power and detuning. It is, however, still very 
large compared to typical probe durations (0.5-2 ms). 

The lifetime of the transverse spin components, however, turns out to be much 
more critical to our experiments. As can be seen from Eq. (}4"9~|) the transverse 
spin components are determined by coherences between magnetic sublevels in 
the x-basis. Therefore anything that affects the lifetime of J x will also affect T 2 . 
In addition, however, the total transverse spin components are also sensitive 
to random phase changes in each atom. As discussed in Sec. 14.31 we can use 
the widths obtained in MORS signals in different experimental settings to 
quantify and separate the effect of different decoherence mechanisms. The 
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Figure 8. Linewidth of the 034 coher- 
ence in the MORS signal as a function 
of optical power. In a single trace it 
is impossible to separate the expected 
power broadening from absorption from 
the light induced collisions. 



Figure 9. Linewidth of the 034 coher- 
ence in the MORS signal as a function 
of atomic density for different optical 
powers. The fact that the slopes are not 
equal reveals light induced collisions. 



FWHM obtained from such signals are related to T2 by: 

r trans [Hz] = — I- (59) 

7Ti 2 [sJ 

We can separate the mechanisms into two main categories: some are medi- 
ated by the probe light and the rest are independent of the presence of the 
probe. Starting with the latter, these combine to a decay rate, Tdark, and con- 
sist mainly of phase changing and spin-flip collisions with the walls and other 
atoms and random phase changes because of atomic motion through inhomo- 
geneous magnetic fields. The effect of these are reduced by the paraffin coating 
on the inside of the glass cells, the diluteness of the atomic sample, and the 
application of additional dc-magnetic fields to cancel field gradients. To deter- 
mine r^ark we measure the width of the mp = 3 «-> 4 coherence for different 
probe powers and find the residual width in the absence of light. An example 
of this is shown in Fig. |H1 As can be seen we obtain a width of the order of 12 
Hz corresponding to a lifetime of the transverse spin of T^ ark w 27ms. We see 
that this will limit but not destroy all correlations between two subsequent 
pulses in e.g. an entanglement experiment as described in Sec. 13.11 

Turning to the probe induced decoherence mechanism, we have already men- 
tioned absorption and subsequent spontaneous emission in the discussion of 
T\. Adding this effect to the other decoherence mechanism, we would expect 
a total decoherence rate of the general form: 

Tidoai = a + b-n + c • P , (60) 

n is the atomic density, P is the light power, and a, b, and c are coefficients, 
which can be determined experimentally. If T is plotted vs. n we would expect 
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a line with constant slope b and offset determined by the optical power. In 
Fig. El we show measurements of the decoherence rate vs. atomic density for 
different optical powers. The results clearly contradict the simple model of 
Eq. (|B0J). since the slope grows with increasing power. It turns out that the 
experiments fit a model: 

Texp — a + b-n + c-P + d-n-P , (61) 

where the size of the expected pure power broadening term, c, agrees with so- 
lutions of the Maxwell-Bloch equations for the full multi-level atomic system 
in the presence of Doppler broadening. The last term could represent light 
induced collisions, but a clear theoretical understanding of the nature of these 
is still missing. For the experimentally relevant densities and powers this term 
contributes around 30 Hz of broadening and is thus the main source of de- 
coherence. We stress that this is a pure T 2 process since we do not observe 
similar features when investigating the decay of the longitudinal spin. Hence, 
the atoms practically decay towards the fully polarized state, i.e. the coherent 
spin state. This was also assumed implicitly in the inclusion of decoherence for 
the entanglement and quantum mapping protocols in Sec. EH where decoher- 
ence was modeled by an admixture of a vacuum state with the same variance 
as the coherent spin state. 

5.3 ENTANGLEMENT RESULTS 

5.3.1 Conditional Entanglement 

We now turn to the experimental demonstration of entanglement generation. 
First the boundary between the classical and the quantum fluctuations has to 
be established. As discussed in Sec. 15. II this projection noise level is found by 
performing several measurements of Xl as a function of the macroscopic spin 
size and verifying a linear increase of the atomic noise of each measurement. 
This linearity combined with nearly perfect orientation of the sample ensures 
the correct projection noise level. Once this is established we implement a 
probing sequence (see Fig. |2J) in which the initial probing pulse is followed by 
a second one after a short delay. To verify entanglement we need to fulfill the 
criterion ([36)1 . in which case our ability to predict the outcome of the second 
probing of the atomic state conditioned on the result of the first measurement 
exceeds the classical limit. 

For the experimental data we calculate the atomic part of the noise by sub- 
tracting the shot and electronics noise of a single light pulse and then normalize 
to the light shot noise level. In Fig. ^] the resulting atomic noise is shown for 
the first and second pulses (with filled and empty circles respectively). Since 
the measurement is of a QND-type and the variance is calculated based on 
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Figure 10. Atomic noise in units of shot 
noise as a function of the macroscopic 
spin size (measured by DC Faraday ro- 
tation). The dotted curve is a quadratic 
fit to the first pulse variances (dots) and 
the dashed curve is the linear part of 
this. Dash-dotted and full-drawn lines 
are fits to the the optimal weight fac- 
tors, a, (x's) and the conditional vari- 
ances (+'s) respectively. 



Figure 11. Deterministic entanglement 
generation. For a fixed atomic density 
negative feedback of the first pulse mea- 
surement is alternately turned on and 
off, thus switching between conditional 
and unconditional entanglement gener- 
ation. 



10,000 independent repetitions of the pulse sequence, the variance of the first 
and the second pulses should be identical. We make a quadratic + linear fit 
to the first pulses' variances (dotted curve) and from this extract the lin- 
ear part, which represent the coherent spin state level. The slope of this is 
0.176(12), which can be compared to the theoretical values 0.187 and 0.165 
obtained from Eq. (fo7)|) with and without the effect of atomic motion included. 
In this fit we used: T=2.0ms, P=5.0mW, and A/2vr = 825MHz. As discussed 
in Sec. El the results including atomic motion are very preliminary, since this 
is an ongoing research topic. The overall agreement with experimental results 
are, however, very encouraging. Note that the small quadratic component is 
caused by various classical noise sources. 

Next we calculate the correlations between the first and the second pulse mea- 
surements as discussed in Sec. 13.11 In Fig. EH the pulses show the conditional 
variance Var(A 2 |A 1 ) + Vea(B 2 \Bi) normalized to shot noise and with shot 
and electronics noise subtracted. The points below the straight line (}3T)j) indi- 
cate that we have in fact created an entangled state between the two atomic 
samples. For the highest densities the noise reduction is up to 30%. The cor- 
responding a-parameters from the minimization procedure ()34j) are plotted 
in Fig. [TUl with crosses. Ideally we would expect these to follow Eq. (J3~9"j) but 
since the atoms decohere the appropriate expression is instead Eq. (|JTJ). Fit- 
ting, we obtain (3 = 0.619(11), which is inserted into Eq. (jlUj) (solid curve). 
Considering the fact that the solid and dash-dotted curves are obtained from a 
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single free parameter, (3, the agreement between experimental results and the 
simple model of decoherence must be considered very satisfactory. In addition, 
we would like to stress the very important point that the two atomic samples, 
which are entangled, are in co mpletely separate environments about 0.5 meters 



apart ([Sherson et all l2005a[ ). This represents a major breakthrough towards 
the creation of truly distant entanglement, which combined with quantum 
teleportation will enable quantum communication over long distances. 



5.3.2 Unconditional Entanglement 



As we have just seen, the results of two probes of the spin state yield corre- 
lated results. The actual results, however, vary from shot to shot, representing 
random realizations of the probability distribution of the spins. That is, we 
create a non-local state with reduced variance but with a non-deterministic 
mean value. Thus, the entanglement only appears when the knowledge gained 
in the first pulse is applied. To create a deterministically entangled state in 
which no knowledge of measurement results is necessary would of course con- 
stitute a very important advance. We have realized this experimentally by 
simply feeding the result of the first measurement pulse back to the atoms 
using an RF-magnetic pulse as discussed in Sec. 14.3.21 This procedure is very 
clo sely related to the wa y in which unconditional spin squeezing is generated 



in (|Geremia et all 120041 ) except that there the feedback is applied continu- 



ously in time, which is more robust against errors in the feedback strength. 



In the experiment we first choose a certain atomic density. Using the linear fit 
of Fig. El we find the projection noise level relative to which entanglement is 
estimated. Next, we vary the strength or the phase of the RF feedback pulse, 
and observe the variance of the second pulse fluctuations as a function of these 
feedback parameters. The feedback giving the lowest variance corresponds to 
the feedback with the optimal gain. In Fig. ^2 we alternate between having 
feedback and not having it. As can be seen clearly, with feedback we obtain 
the same variance for the second pulse as we do without feedback for the con- 
ditional variance. This means that we have in fact created a state, which is 
deterministically entangled. That is, the state after the first probe and sub- 
sequent feedback of the measurement result has zero mean and a variance 
reduced by 19% compared to the projection noise level. When feedback is 
turned on the variance of the second pulse and the conditional variance coin- 
cide, which means that knowledge of the measurement result cannot improve 
the degree of entanglement beyond the entanglement proved by the second 
pulse measurement alone. 
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5.4 QUANTUM MEMORY RESULTS 



We will now present the ex perimental demonstratio n of quantum memory, 
which has been published in f Julseaard et all |2004a ). The used protocol was 



discussed in Sec. 13.21 and, briefly, involves 1) preparation of the initial atomic 
state in a coherent state via optical pumping, 2) mapping of one of the light 
quadratures through the off-resonant Faraday interaction, Eqs. I30al -d. and 3) 
mapping of the second light quadrature by a direct measurement and subse- 
quent feedback (an illustration of the timing sequence is shown in the inset 
of Fig. H2J). For the reasons discussed in Sec. 13.2.31 we have not been able 
to retrieve the mapped state. Instead we have performed a destructive re- 
construction of the mapped state. This is done by waiting for a time r and 
then sending a readout light pulse through the atomic sample. Measuring 
the X-component of the outgoing light will then give information about the 
stored atomic P-component according to Eq. (|30a|) . If a 7r/2 rotation in the 
atomic XP space is performed prior to the readout pulse we obtain informa- 
tion about the stored atomic X-component. Repeating this 10,000 times, the 
statistics for the atomic variables after the storage procedure can be recon- 
structed. The first thing to check is that the mean value of the stored state 
depends linearly on the mean value of the input light state. This is shown 
in Fig. First of all we note that the linear dependence is clear for both 
quadratures. This completes the proof of classical memory performance. The 
next thing to note is that the slope is not unity, which means that the stored 
state has a different mean value than the input state. The reason for this is 
discussed further below. For the quadrature mapped straight from the back 
action of the light onto the atoms we have the gain of #ba = 0.836 and for the 
quadrature mapped via the feedback we have g-p = 0.797. In order to verify 
quantum storage we also need to consider the shot-to-shot fluctuations in the 
stored state, which for a Gaussian state are fully characterized by the vari- 
ance of the state. The experimentally reconstructed variances of the atomic 
quadratures, that is the variance of the readout pulse with the one unit of 
shot noise intrinsic to the read-out pulse subtracted, are shown in Fig. ED 
Also shown is the ideal quantum limit on the variance for a perfect mapping 
and the classical limit for uq = 4. As can be seen the variance is more or less 
independent of the mean value of the input light quadratures. The fidelity of 
the stored state for a Gaussian distribution of input states with mean photon 
number no can be calculated given the measured gains and variances, a x and 
0"p, according to: 

(62) 



(2n (l - #ba) 2 + 1 + a x )(2n (l - g F ) 2 + 1 + a P ) 

With the experimentally measured values we get F = (66.7± 1.7) % for n = 4 
and F = (70.0 ± 2.0)% for n = 2. The best classical fidelity was recently 
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Figure 12. Mean value of the read out 
pulse as a function of the mean values of 
the input light variables, X™ and P™, to 
be stored. Inset: the strong classical and 
weak quantum pulses in opposite polar- 
izations. Between the input and the out- 
put pulses are the feedback pulse and 
the optional ir/2 pulse. 
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variances for ran = 4 



derived by Hammerer et al. ( 2005b! ) 
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for coherent states drawn out of a Gaussian distribution with mean photon 
number n. This means that -F c i ass decreases monotonically from unity for the 
vacuum state to 1/2 for an arbitrary coherent state. For n = 4, 2 we obtain 
the classical boundaries of 55.6% and 60.0% respectively. This verifies that the 
storage of the light state in fact constitutes a quantum mapping. The results 
shown were obtained for a pulse duration of 1 ms and the memory has been 
shown to work for up to 4 ms delay between the two probe pulses. Note that 
we have chosen to calculate the fidelity as the average of the squared overlap 
between the stored state and the ideally stored state. For non-unity gain this 
decreases very rapidly with coherent states having large amplitudes. However, 
one could argue that a storage with an arbitrary but known gain constitutes 
just as useful a memory as unity gain memory. If analyzed solely in terms of 
the added noise, our memory would therefore perform much better than the 
previously stated results, which can therefore be viewed as a lower bound on 
the memory capability. 



Note also that the choice of the optimal gain which maximizes the fidelity 
depends on the class of available states. For classes of coherent states limited 
in their amplitudes discussed in this section, the experimental gains quoted 
above are, in fact, close to optimal ones. 
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Measured feedback gain 

Figure 14. Reconstructed variances of X^ ut (o's) and P£ ut (x's) as a function of feed- 
back gain. The curves are theoretical curves based on independently determined 
noise parameters. 

5.4.1 Decoherence 

In Fig. El we show the final atomic variances as a function of the feedback 
gain. These are compared to the variances expected from Eqs. (I44|) with in- 
dependently determined decoherence values of ft = 0.61 and ( = 0.75. We 
stress that this is not a fit to the data, which means that we understand the 
level of experimentally determined atomic noise quite well. The figure clearly 
shows that, because of the decoherence and light loss, if the feedback gain is 
increased towards unity the noise grows dramatically. With this, the fidelity 
quantified by Eq. (jS2j) can be optimized with respect to the feedback gain. As 
can be seen from the values of ft and £ the light loss and atomic decoherence 
are significant. The high light loss is due to the fact that the glass cells contain- 
ing the atomic vapor were not anti-reflection coated. Therefore each glass-air 
interface contributes about 4% loss. The main source of atomic decoherence 
is light assisted collisions, which change the phase of the atomic coherence 
without affecting the macroscopic spin size J x . The atoms are driven towards 
a coherent state, which justifies the use of the simple model of beam splitter 
admixture of vacuum. 



6 Conclusions 



We have described the progress in development of the deterministic quantum 
interface between light and macroscopic atomic ensembles. Theory of the inter- 
face is based on canonical variables which provide a convenient language when 
different physical systems, such as light and atoms, are considered. Surpris- 
ingly enough, simple dipole off-resonant interaction of light with spin polarized 
atomic samples with high optical density provides a powerful tool for quantum 
state engineering and transfer. To perform sophisticated quantum information 
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protocols, the interaction can be combined with a quantum measurement on 
light and subsequent feedback to the atoms. Two central experiments of this 
review, entanglement of two macroscopic objects and quantum memory for 
light, are performed following these general steps. 

Future perspectives include multipass protocols towards more efficient quan- 
tum memory including the retrieval process, as well as various types of quan- 
tum teleportation involving atomic ensembles. Another exciting future direc- 
tion is demonstration of deterministic qubit state manipulation with the tool 
box described in this paper. It is appropriate to stress that the language of 
canonical variables is fully applicable not only to Gaussian states but to all 
single mode states, including qubit and Fock states. This language is asso- 
ciated with homodyne measurements in quantum optics. Addition of single 
photon counting to the experiments described in this article may pave the 
road towards more efficient and robust quantum information processing and 
communication. 
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A Effect of Atomic Motion 



Since in our experiments the atoms are at room temperature and, for exper- 
imental reasons, the light beam does not cover the entire cross section of the 
atomic sample, the atoms move across the beam several times (~ 10) during 
the time of a pulse. This averaging effect insures that all atoms spend roughly 
the same amount of time inside the beam but, as we shall see, it still has 
important implications for the noise properties. In brief, the atomic motion 
modifies the projection noise level and acts as an additional source of deco- 
herence since two subsequent probe p ulses interact with the atoms d ifferently. 
The results are related to the work of iKuzmich and Kennedy! ((2004) . 
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A.l MODELING ATOMIC MOTION 



To be more quantitative, we introduce new pseudo-angular momentum oper- 
ators J q — > J2iPiFq l \ where Pi is the fraction of time the i'th atom spends 
interacting with the laser beam and q — x,y, z. These have the commutator: 



l-Wat JVa 



N^t ATa- 



Eft^E^ = Eft 2 W .^! ] = <E^i i} - (a- 1 ) 

_i=l i=l J i=l i=l 

This leads to the Heisenberg uncertainty relation (for a highly polarized sample 
with F x « F) 

/A^at \ /A/at \ / T x 2 

Var (Eft^J Var (Eft^J > (^ 2 (! + ff2 )) ( A -2) 

where we have introduced the mean p = (p^ and variance Var(pj) = a 2 -p 2 
of With this definition, a is the relative standard deviation of p. Since, for 
the coherent spin state, 

/Af at \ T 

Var f E^i ) = 2 p2(1 + = Var ( CSS )> ( A - 3 ) 

this highly polarized state corresponds to a minimum uncertainty state. The 
measured noise is then limited by the Heisenberg uncertainty principle and we 
confidently call this projection noise. To maintain the correct commutation 
relation [X, P] = i we experimentally normalize the atomic operators to the 
measured projection noise, i.e. instead of defining X = J y j\fj we effectively 

have X = ZiPiFW/y/W + a 2 ). 

The average fraction of time p each atom spends inside the beam is clearly 
p = -Abeam/Axsii- Let us now discuss the scaling of a 2 with simple physical 
parameters. The fact that the variance may be non-zero arises from the finite 
time available for the averaging process carried out by the atomic motion. A 
typical traversing time across the vapor cell is r = L/v where L is the cell 
dimension and v is e.g. the one-dimensional rms speed of the atoms. We may 
think of this atomic motion as n independent journeys across the vapor cell 
volume, where n ~ T/r = Tv /L. We then model the motion through the 
beam with mean occupancy p by assuming in each walk across the cell volume 
that either (1) the atom spends all the time r inside the beam, which should 
happen with probability p, or (2) the atom spends all the time r outside the 
beam which should happen with probability 1 — p. 

We then count the number of times n; ns id e that an atom was inside the beam 
out of the possible n journeys. In this simple model n insi d e is a stochastic 
variable which is binomially distributed with mean np and variance np(l —p). 
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We are interested in the fraction of time (~ n- ms ide/n) spent inside the beam. 
It follows (n insidc /n) = p and a 2 = Va.r([ni nsi( ±e/n}/p) = (1 — p)/np. Hence the 
simple model leads to 

^4bcam -, 2 (-^cell ^.bcam)-^ /. .\ 

p = — and a = — (A.4) 

^icell ^ibcam-f "^0 

Note the characteristic scaling with T _1 and with the area (A ce n — A heam ) not 
covered by the light beam (when A^, eSuia is close to its maximum value A ce n). We 
note that due to the simplicity of the above model the absolute numbers should 
only hold as an order of magnitude estimate. Numerical simulations performed 
for a cubic cell have shown that the relative variance a 2 is roughly four times 
smaller than the estimate above. Also, due to the Doppler broadening, the 
effective detuning differs from atom to atom and causes an increase in a 2 . 



A.2 ATOMIC MOTION AS A SOURCE OF DECOHERENCE 



To see how atomic motion acts as an effective source of decoherence, imagine 
that we perform some manipulations of atoms using one laser pulse and subse- 
quently probe these manipulations with another laser pulse. Since atoms move 
during interactions the probed quantum operator changes in time. Comparing 
the operator at the 1st and 2nd times we get 

( N&t c\ N&t f\ \ Nllt 

Var ^Pi,2nd^j,2n d ~ Pi,^ F zlst = Va K F i,lst) ((Pi,2nd - Pi,lstf 
\i=l i=l / i=l 

= i . 2p 2 a 2 = 2Var(CSS)(l - (3) with [3 ' 



1 + a 2 

(A.5) 

We assumed Pi t i s t and Pi^nd to be uncorrelated, which is reasonable since a 
collision with the cell wall randomizes the velocity of the atoms. Also, we took 
^i'lst = -^2,2nd- This corresponds to having no decoherence at all apart from 
the effect of atomic motion which is the only effect studied in this calculation. 
To interpret the above calculations we consider a standard decoherence calcu- 
lation. Consider a true spin operator J z subject to decoherence parametrized 
by the number (3 such that 

J g ,iBt -»■ Jz,2nd = PJ z ,ist + yfi- /9 2 J vac with Var(J vac ) = J - = Var(CSS). 

(A.6) 

Then the operator changes by an amount characterized by the variance 



Var (jf d - J l z st ) = Var (j 2 lst (l - (3) - y/l - (3 2 J vac ) 
= J(l -/3) = 2Var(CSS)(l - (3) 



(A.7) 
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which is exactly the same as in (|A.5[) . We are led to the conclusion that 
atomic motion inevitably gives rise to an effective decoherence. We thus see 
that, whereas the increased coherent spin state noise with increased a 2 might 
seem to suggest that non-classical states are more easily created (by producing 
states with noise lower than the CSS), this is compensated for by an increased 
decoherence of the state. Therefore, higher a 2 does not lead to higher fidelity 
protocols. 



B Technical Details 

B.l LIGHT POLARIZATION AND STARK SHIFTS 

Let us calculate the Stark effect from the probe laser on the magnetic sublevels 
\F,m). We let the light be strong and linearly polarized along the vector 

ei = e x cos a + e y sin a, (B.l) 

i.e. a is the angle between the macroscopic spin direction (the x-axis) and 
the probe polarization direction. Light is propagating in the z-direction. The 
Stark effect on magnetic sub-levels is much weaker than the splitting caused 
by the constant bias magnetic field and can be calculated in non-degenerate 
perturbation theory from the interaction Hamiltonian Q. The a® term is 
common to all levels, the a± term is zero on average since (Sz) = 0, and we 
are left with the higher order components proportional to 0,2 ■ 

With di being the creation operator along the strong direction we have d x = 
ax cos a and d y = d\ sin a (neglecting the orthogonal direction to d\ which will 
be in the vacuum state for linear polarization). With S± = S x ±iS y we derive 
from (gj) that 

(S + W) = ^e- and (<L( f) ) = M e -- (B.2) 

where 4>(t) is the photon flux and Stokes operators are normalized to photons 
per second. In order to calculate the higher order terms of the interaction 
Hamiltonian for a single atom we leave out the integral / . . . pAdz in (jZJ) and 
renormalize light operators (by absorbing the speed of light c) to photons per 
second and find 

**** = ~8lA2^ a2 ' m ' " [fx ~ $ C ° s(2a) " [jx3y + jy3x] Sin(2a) ) • 

(B.3) 
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We also replaced j± by j x ± ij y . We need to calculate the expectation value of 
this Hamiltonian for the different energy eigenstates \m) quantized along the 
x-direction. We have 



(m\ jI \m 



(m\ f \m 



(m\ j 2 z \m 



(m\ jxjy + jyjx \m 



m 2 , 



F(F + 1) -m 



2 



F(F+1) - m 



0. 



2 2 (B.4) 



The first of these is obvious since \m) is quantized along the x-axis. We have 
(pv) = symmetry and the value is found from the fact that j 2 , + 

Jy + Jz = F(F + 1). Also, by symmetry we have in an eigenstate of j x that 
{m\j y j x \m) = m ■ (m\ j y \m) = which leads to the final line. Calculating the 
expectation value of (jB.3|) we get 

l + 3cos(2a;) 2 1 + cos(2a) . 
— - • m z — L F(F + 1) 



(B.5) 

What is important for us is the shift Su m+lt7n = {E^^ — E^ &rk )/h of a 
magnetic resonance line which then becomes 



^SJHz] = • • [1 + 3 cos(2a)] • [2m + 1]. (B.6) 



This Stark shift is problematic for several protocols, especially the setup with 
two oppositely oriented samples. Note, that for atoms polarized in the mp = 4 
state the relevant transition = 4 <-> mp = 3 has a Stark shift proportional 
to 2 ■ 3 + 1 = 7. An oppositely oriented sample with mj? = —4 has for the 
transition mp = —3 <-> = —4 a Stark shift proportional to 2 ■ (—4) + 1 = 
—7. Hence, these two transition frequencies cannot be overlapped both in the 
presence and absence of light (see Fig. IB.l| for an illustration). 

One remedy for this is to choose the light to be linearly polarized at an angle 
a = 54.7° corresponding to 3cos(2a) = —1, and the Stark term disappears. 
As we shall see in Sec. IB.2l this gives rise to other problems. 

Another remedy is to add an extra bias magnetic field along the x-direction 
when the laser light is on. In this way the frequency of the desired transitions 
can be kept stable. This is the approach we have taken and it works well. One 
should note though, that with our laser pulse timing in the vicinity of one 
millisecond, it is not completely trivial to make a magnetic pulse following 
the laser intensity since eddy currents in metallic parts and induced electric 
fields in other current loops for magnetic fields should be taken into account. 
A convenient diagnostic method is to apply a classical shift to the spin states 
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With light 
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Figure B.l. An illustration of the problem with the light induced Stark shift of the 
Zeeman sublevels. Without applying additional fields the two important Zeeman 
resonances cannot be overlapped both without light (a) and with light (b). 

along the lines of Eq. (jlHj) prior to the application of two laser pulses. An S y 
detection of the light will show in real time the mean value of the spin state 
components J' y and J' z . These should be constant through each laser pulse and 
conserved in the dark time in between (apart from a possible decay) if all 
frequencies are well adjusted. 



B.2 INFLUENCE OF LASER NOISE 



In this section we make a brief comment on laser noise entering the atomic 
samples. This is by no means a complete analysis but it describes some im- 
portant issues in connection with the choice of laser polarization. We start 
out with the interaction Hamiltonian (J7J) and also assume that the atoms are 
polarized along the x-axis and placed in a bias magnetic field, leading to the 
magnetic Hamiltonian H mag = Qj x . This leads to the following equations of 
motion for the transverse spin components j y and j z : 

- h (z, t) = -n Jz + — - {-a lSz3x ^ (B.7) 

+ a 2 (-(2§ x + $) [j x j z + j z j x ] - 2S y [j z jy + 3yj z })} , 

-j z (z, t) = Qj y + (l-V, J,./,, + jyj x ] - AS y \jl - f y )} . (B.8) 



The Larmor precession terms take all interesting dynamics to the frequency 
component around Q. Let us see which terms above couple ^-components 
of light into the spin state. First, the a,iS z ] x term is our favorite interaction 
term (JT5j) used in all quantum information protocols. It consists of an atomic 
operator j x which is constant equal to ±F for all practical purposes. This is 
multiplied by S z , whose il-components drive dj y /dt. 
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In Eq. f)B.8|) the final term is proportional to j 2 — j 2 . Practically, this is also 
constant equal to F 2 — F/2. It is then multiplied by S y , whose f2-components 
drive dj z /dt. This is an unwanted effect. The ratio of the unwanted to wanted 
noise is found by squaring these contributions. We get 

Unwanted noise p 2 /a 2 \ 2 Noise 



Wanted noise \aiJ Noise (S 

For our typical values of detuning we have 02/01 ~ 10~ 2 and the above ratio 
becomes ~ 0.02 ■ Noise(S , ?/ )/Noise(5'^) for F = 4. If our laser beam is polarized 
along the x- or ?/-axis with a clean linear polarization, the noise of S y and S z 
at frequency Q can be shot noise limited, i.e. by quantum noise (amplitude 
noise of the laser does drive into S y and S z in the case of clean linear polariza- 
tion). In this case the unwanted noise only contributes ~ 2% of the total noise 
pile up. But if we choose arbitrary polarization directions in the xy-plane, the 
S^-component will have non-zero mean value, and the fluctuations at Q will 
essentially be the amplitude noise of the laser at Q. In this case, to keep the 
last term of Eq. (jB.8|) from accumulating extra noise, one requires the laser 
intensity to be shot noise limited at Q (which is a more difficult condition to 
meet than clean linear polarization). We thus have one motivation for choosing 
the laser to be polarized along the x- or ^-direction and not in between. We 
typically do this in our experiments and for this reason we have to compen- 
sate the Stark splitting discussed in Sec. IB.ll For a further d i scussi on of the 



different higher order terms in the interaction, see (|JulsgaardL 120 03). A more 
thoroug h discussion of quantum n oise with the higher order terms included is 
given bv IKuprivanov et al. (|2005h . 
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